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We use rigorous group-theoretic techniques and molecular dynamics to investigate the connection 
between structural symmetry and ionic conductivity in the garnet family of solid Li-ion electrolytes. 
We identify new ordered phases and order-disorder phase transitions that are relevant for conduc¬ 
tivity optimization. Ionic transport in this materials family is controlled by the frustration of the Li 
sublattice caused by incommensurability with the host structure at non-integer Li concentrations, 
while ordered phases explain regions of sharply lower conductivity. Disorder is therefore predicted 
to be optimal for ionic transport in this and other conductor families with strong Li interaction. 
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Solid-state Li-ion batteries have generated keen inter¬ 
est as advanced energy storage systems with applica¬ 
tions in portable electronics and electric vehicles [1]. A 
highly conductive and stable solid ceramic ionic conduc¬ 
tor can enable high-energy density by serving as a pro¬ 
tection layer for Li-metal anodes, and a stable electrolyte 
for high-voltage cathodes. At the same time, replacing 
flammable organic electrolytes with an inorganic com¬ 
pound will remove the threat of thermal runaway catas¬ 
trophes. In this context, solid-state inorganic Li con¬ 
ducting oxides with the garnet structure are currently 
considered to be among the top promising candidates for 
solid electrolytes due to a combination of their perfor¬ 
mance and stability properties. The garnet-like conduc¬ 
tor Li 5 La 3 M 2 0 12 (M = Nb, Ta) was first reported by 
Thangadurai et ai [2] and to date an ionic conductivity 
as high as a = 3 x 10 _3 5'/cm has been measured [3] for 
the cubic Li 7 La 3 Zr 2 0 12 garnet. The prototypical garnet 
structure is cubic (Ia-3d No. 230) with the chemical for¬ 
mula LUb 3 C 2 0 12 . This structure is remarkably robust in 
regard to changes in cation composition (B=La, Ca, Ba, 
Sr, Y, Pr, Nd, Sm-Lu and C=Zr, Ta, Nb, Nd, Te, W) and 
is able to accommodate Li concentrations x ranging from 
3 to 7. There are two distinct crystallographic cation sites 
that Li can occupy: octahedral (oct) 4^9 and tetrahedral 
(tet) 24d. Other cations occupy the 24c (B) and 16a (C) 
sites, and oxygen occupies the 96h site. It has been rec¬ 
ognized that Li site arrangements vary with Li content 
x , but there are conflicting reports of the trends coming 


from experimental measurements [4, 5]. The complex¬ 
ity presents a difficulty in understanding both structure 
and ionic transport mechanisms. Unfortunately, X-ray 
diffraction can provide limited help in characterizing Li 
arrangements due to the weak scattering from Li, and 
one has to resort to complicated neutron scattering meth¬ 
ods. The nature of ionic conduction is even less clear, as 
it is nearly impossible to observe directly, while at the 
same time a complex relationship between crystal struc¬ 
ture and conductivity is evident. Multiple phases have 
been reported across the composition space with widely 
ranging conductivity values, which depend not only on 
stoichiometry but also preparation routes. For example, 
it was reported that Li 7 La 3 Zr 2 0 12 is tetragonal (t-LLZ) 
with two orders of magnitude lower conductivity than 
the cubic (c-LLZ) polymorph [6]. Despite the promise 
and strong interest in this material there is limited un¬ 
derstanding and consensus regarding the structure and 
transport mechanisms, particularly regarding the opti¬ 
mal Li vacancy concentration and site occupation [7]. 

In this work, we present a novel analysis method and 
use it to derive a comprehensive description of the atomic 
mechanisms governing garnet crystals, which is able to 
explain, in a predictive way, both structural and dy¬ 
namic effects in the wide composition space. Our uni¬ 
fying description of this important materials class is en¬ 
abled by a novel combination of computational methods 
that include systematic hierarchical space group sym¬ 
metry analysis coupled with high-throughput atomistic 
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computations and molecular dynamics. Structural opti¬ 
mization and temperature dependent transport was in¬ 
vestigated with both Born-Oppenheimer ab-initio molec¬ 
ular dynamics (AIMD) using Quantum Espresso [8] and 
classical molecular dynamics (CMD) using LAMMPS [9]. 
Details of the computational procedures are given in the 
Supplementary Information. We note that our analy¬ 
sis focuses only on single-crystal transport and neglects 
effects of grain boundaries which can be significant in 
sample measurements. 

The arrangement of Li ions on the available crystal¬ 
lographic sites is arguably the main degree of freedom 
that controls both structure and transport characteris¬ 
tics of the garnet solid electrolytes. In order to identify 
the ground state Li orderings at each composition, we 
need to compute energies of all possible orderings. How¬ 
ever, the complexity of even the primitive unit crystal 
cell makes it impossible: in a simple cubic representation 
of LLZ there are 72 Li sites, of which 56 are occupied, 
resulting in 10 20 possible arrangements. In order to deal 
with the enormous combinatorial space of configurations, 
we combine rigorous group-theoretic analysis for classi¬ 
fying Li site arrangements for each composition together 
with automated computations of total energies, in order 
to identify occupancy rules. The resulting complete map¬ 
ping of lowest energy structures at each composition is 
then used to predict new ordered or disordered phases 
of garnets and to postulate the dependency of conduc¬ 
tivity on composition. The first step in our analysis is 
the determination that the occupancy ratio of Li among 
tet and oct sites is universal: it depends only on the Li 
concentration x in the composition. To demonstrate this, 
for each value of x we generate a set of ^100 candidate 
configurations of various oct:tet occupancy ratios. Each 
structure is relaxed using DFT and the occupancy of the 
lowest energy configuration is automatically extracted 
using a topological partitioning algorithm. This proce¬ 
dure involves partitioning of the crystal by constructing 
a compact geometric network of non-overlapping polyhe- 
dra with O atoms as the vertices. Each polyhedron is 
then analyzed at each time step by a 3D convex hull al¬ 
gorithm to determine whether there is a Li ion present 
within the boundaries. The advantage of this topology- 
based approach is that it gives accurate indications of site 
occupancies, free from thermal noise, without the need 
of arbitrary distance-based metrics. 

The obtained occupancy trend is surprisingly linear, 
with Li favoring tet site occupancy when the Li con¬ 
tent is low (x < 4.5). While qualitatively similar 
trends have been observed in several experimental re¬ 
ports [4, 5, 10, 11] quantitative results show a significant 
variation caused by experimental uncertainties and differ¬ 
ences in material preparation. In order to investigate the 
effects of dynamics on this trend we performed AIMD 
at several temperatures. Fig.l shows the temperature 
dependent occupancy over different garnet compositions 


averaged over each AIMD run. It can be inferred that 
the trend is largely independent of temperature. The 
weak temperature dependence is in agreement with sev¬ 
eral in situ neutron diffraction experiments conducted 
on Li 5 La 3 Ta 2 0 12 [12] and Li 6 BaLa 2 Ta 2 0 12 [13] between 
room temperature and up to 900 K that reported negli¬ 
gible variation in the lithium distribution across the tet 
and oct sites. 
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FIG. 1. Linear trend of Li occupancy (black) of oct and 
tet sites as a function of Li content. Temperature de¬ 
pendent oct (green) and tet (red) Li occupancy (500-1300 
K) do not significantly deviate from occupancies at 0 K. 
Labels denote: c - cubic, LLZ - Li 7 La 3 Zr 2 0 12 , LLZA1 
- Li6.25Alo.75La 3 Zr 2 0 12 , LLZT - Li 6 La 3 ZrTa0 12 , LLT - 
Li 5 La 3 Ta 2 0 12 , LNT - Li 3 Nb 3 Te 2 0 12 

The occupancy constraint allows us to reduce the num¬ 
ber of candidate configurations in the search for the 
ground state at each composition (3 < x < 7), but the 
number of possibilities is still on the order of 10 19 . To 
reduce the number of relevant configurations we perform 
a systematic search through crystallographic symmetry 
groups, subject to the universal constraint of the site oc¬ 
cupancy ratio. This mathematical procedure enables us 
to efficiently analyze the entire space of possible ordered 
ground state structures in the garnet crystal family, using 
only a few hundred configurations. The group theoretic 
analysis procedure starts with crystallographic projec¬ 
tions of the prototypical cubic symmetry group Ia3d (no. 
230) of the host lattice onto lower-symmetry subgroups 
in order to enumerate ordered structures at partial Li oc¬ 
cupancies. The subgroup projection is implemented us¬ 
ing tools provided by the Bilbao Crystallographic Server 
[14]. We only consider cubic and tetragonal candidate 
subgroups, motivated by the experimental observation of 
only these two lattice types in the garnet family. First, 
for each Li concentration x, we determine all subgroups of 
group No.230 with the k-index up to 2. (The k-index indi¬ 
cates the multiplication factor relating the volume of the 
primitive cell of the subgroup with respect to the primi¬ 
tive cell of the original prototype structure.) This limit is 
imposed in order to limit the search to ordered structures 
with a reasonably compact description of a highly ordered 
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phase. We then analyze the splittings [15] of Wyckoff po¬ 
sitions in each subgroup and conjugacy class in order to 
identify symmetrically related sets of Li sites and their 
multiplicities. The oct: tet occupancy ratio constraint is 
introduced by requiring that each resulting structure ac¬ 
commodates the correct ratio with the available Li site 
multiplicities. Finally, for each ordered structure that 
satisfies these constraints we perform a DFT variable-cell 
optimization calculations and identify the Li configura¬ 
tion of lowest energy for each composition x. The main 
result of this analysis is the non-trivial finding that for 
each of the integer values ofx = 3,4,5,6,7 there exist or¬ 
dered configurations with cubic or tetragonal symmetry, 
and that some are more robust than others with respect 
to the tendency to disorder. 



FIG. 2. 2D representation of the ground-state Li arrangement 
for x as 3 (left), x — 6 (middle) and x — 7 (right). Tet sites 
are denoted by grid intersections, and oct sites occupy the 
segment midpoints. The other constituents of the crystal are 
omitted for ease of viewing. 

The topology of the Li site network determines the 
degree of ordering and consequently the transport char¬ 
acteristics in a non-trivial manner. Each tet site shares 
a face with four neighboring oct sites while each oct site 
shares a face with two neighboring tet sites. In order to 
easily visualize the Li arrangements, we map the topology 
of the 3D Li network onto a 2D grid (which provides an 
only locally correct mapping). Fig. 2 shows the 2D rep¬ 
resentations of the resulting ground-state configurations 
of Li concentrations x=3, 6 and 7. In the following we 
discuss in detail the Li configurations we identify for each 
integer concentration x. We note that it is not possible to 
find ordered structures (within the specified constraints) 
for off-integer concentrations, so they are predicted to be 
cubic and disordered. 

x=7. From the subgroup projection procedure, the 
highest symmetry subgroup that allows for the occu¬ 
pancy ratio of 1:6 (tetioct) is the tetragonal group 
I4i/acd (No. 142). It is not possible to find any cubic 
subgroups of group No. 230, subject to the constraints 
of occupancy and unit cell size. DFT calculations iden¬ 
tify the tetragonal layered structure as the lowest energy 
configuration. In this configuration Li occupies all the 
oct sites as well as 1/3 of the tet sites in a characteristic 
well-ordered layered arrangement (illustrated in Fig. 2 
right), which imparts a noticeable tetragonal distortion 
(c/a=1.04 for t-LLZ) to the unit cell. XRD/Neutron 


diffraction results [6] confirm the tetragonal nature of 
the ground state structure of LLZ, as well as the lay¬ 
ered Li atom arrangement. Li in an oct site (489 in the 
prototype structure), is distorted to occupy the Wyckoff 
position 96h when one of the neighboring tet sites is oc¬ 
cupied. This oct site distortion is present in other garnet 
compositions as well, indicating a significant degree of 
interaction among Li in the sublattice. The tetragonal 
structure is quite stationary, which is indicated by the 
high energy cost to move a Li ion from an oct site to 
a tet site (^0.7eV) [16] and the tendency to rearrange 
subsequently back to the original state. In addition to 
the tetragonal structure, however, we find multiple low- 
symmetry configurations with the Li occupancy of 1.5:5.5 
(tet:oct), which deviates slightly from the expected oc¬ 
cupancy trend. These structures are higher in energy 
by as little as 5 meV per Li site and the cell shape is 
much closer to cubic, as expected from lower symmetry 
configurations. This is consistent with experimental ob¬ 
servations that a phase transition can turn the ordered 
tetragonal LLZ (t-LLZ) to a disordered cubic form (c- 
LLZ) [17]. 

x=6. The occupancy trend predicts a tet:oct ra¬ 
tio 1.5:4.5 for this composition. A representative com¬ 
pound of this composition is Li 6 La 3 ZrM0 12 , where M = 
Nb 5+ (LLZN) or Ta 5+ (LLZT). Symmetry subgroup anal¬ 
ysis of Wyckoff position splittings identifies only one pos¬ 
sible cubic group No. 198 (P2i3). Within this space 
group the half-half mixing of cations (Zr:Ta / Zr:Nb) can 
also be accommodated in several ways. Overall there 
are 48 structures corresponding to arrangements of Li 
ions and cations within this space group, and we per¬ 
form DFT variable-cell relaxation calculations to identify 
the ground state configuration for both mixed Zr-Ta and 
Zr-Nb cation compositions. The lowest energy structure 
exhibits a rock-salt type ordering of the 4+/5+ cations, 
and a simple pattern of Li ions, where each tet Li has 
exactly 3 oct neighbors (see Supporting Information). In 
addition, each tet Li has zero tet neighbors, and exactly 
half of the tetrahedral sites are occupied, (see Fig. 2 
middle). The excitation energy to disturb the ordering is 
calculated to be 23 meV/Li site. Although experimental 
results report that x = 6 compositions are indeed cubic 
[18], a detailed characterization of Li ion arrangement is 
presently lacking. Our finding presents new information 
explaining the likely ground state structure of this com¬ 
position, which can anticipate confirmation from neutron 
diffraction studies of carefully annealed samples. 

x=5. The prescribed occupancy ratio of 2:3 (tet:oct) 
is possible to obtain within the tetragonal group No. 98 
(I4]_22) as the highest symmetry compatible with the Li 
site Wyckoff position splitting. Group No. 80 (I4i) is a 
tetragonal subgroup of No. 98, and also accommodates 
the required occupancy. DFT variable-cell relaxation cal¬ 
culations of Li 5 La 3 M 2 0 12 (M = Nb 5+ or Ta 5+ ) identify a 
tetragonal structure of group No. 98 as the lowest energy 
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configuration. However, we also find multiple tetragonal 
structures corresponding to space group No. 80 that are 
only 9 meV per Li site higher in energy than the ground 
state. We predict that this degeneracy is likely to yield 
a statistically averaged cubic structure at room tempera¬ 
ture with no long range Li order, as synthesized, for this 
composition [19]. We propose that it may be possible to 
detect tetragonal domains with sensitive neutron diffrac¬ 
tion equipment, for samples annealed at lower tempera¬ 
tures. Similar to the x = 6 case, we find that the rela¬ 
tive energies of arrangements inLi 5 La 3 Ta 2 0 12 (LLT) and 
Li 5 La 3 Nb 2 0 12 (LLN) garnets are very close, predicting 
similar structural and transport behavior in these com¬ 
positions. 

x=4. The tetioct occupancy ratio of 2.5:1.5 can be 
accommodated by Wyckoff position splitting in tetrag¬ 
onal groups No. 117, 116, and 95. By computing total 
energies for each of 142 possible optimized configurations 
we find that the lowest energy structures belong to space 
group No. 95. The one identified as the lowest energy 
ordered structure is detailed in the Supporting Informa¬ 
tion. The next-lowest energy structure in this family is 
higher in energy by 20 meV/Li site, which also makes it 
challenging to observe especially if the high-temperature 
synthesis procedure involves a rapid quenching step. Un¬ 
less carefully annealed, this composition is likely to be 
cubic and disordered, as synthesized. 

x=3. In the classical silicate garnet structure (space 
group No. 230) all the 24d sites (tet) are occupied but 
none of the 4$9 (oct) sites. In the case of Li garnets, our 
subgroup projection approach identifies a unique struc¬ 
ture with group No. 230 as the highest symmetry com¬ 
patible with full tet occupancy (as dictated by the occu¬ 
pancy trend in Fig. 1). This well-ordered cubic structure 
is strongly favored energetically (shown in Fig. 2 left), 
requiring more than 1.0 eV of energy to displace a Li 
ion from its tet position. Previous neutron diffraction 
studies [20] have determined the lithium distribution in 
Li 3 LaTe 2 0 12 and Li 3 Nd 3 W 2 0 12 and conform to our pre¬ 
diction. 

The overall picture is that both tetragonal and cubic 
phases appear in the composition range, but the ionic 
conductivity depends primarily on the degree and ro¬ 
bustness of order in the Li sublattice. Results on Fig. 
3 indicate that the conductivity, as a function of compo¬ 
sition, exhibits a non-trivial trend where it is strongly 
suppressed by the presence of a stable ordered Li ar¬ 
rangement. In the structure analysis above we identi¬ 
fied x = 3,4, 6 and x = 7 (t-LLZ) as having well-ordered 
stable ground state configurations in the Li sub-lattice 
(Fig. 2) with non-negligible excitation energies. Li con¬ 
duction within an ordered sublattice involves creation of 
defects, which costs energy. On the other hand, much 
shallower energy landscapes results in better ionic mo¬ 
bility in disordered configurations that arise in situations 
of frustration. These are exactly the compositions, such 


as half-integral x, where it is mathematically impossible 
to find a simple ground state with a symmetry arrange¬ 
ment compatible with the host crystal lattice. As a re¬ 
sult the Li network is intrinsically frustrated, disordered 
and mobile. We note that the interplay between possi¬ 
ble symmetries of carrier configurations and conductivity 
is reminiscent to the order-disorder controlled transport 
observed in frustrated models, such as the triangular 2D 
anti-ferromagnetic Ising model [21]. 



Li concentration (x) 

FIG. 3. Ionic conductivity (S/cm) of Li + for temperatures 
500 K to 1700 K as a function of Li content x in the garnet 
crystal, from CMD simulations. At x — 7, the conductivity of 
the tetragonal (t-Li7) (denoted as hollow squares) and cubic 
(c-Li7) phases of Li 7 La 3 Zr 2 0 12 have been depicted. 

In the case of x = 6, the newly identified stable or¬ 
dered ground state structure results in poor mobility at 
lower temperatures. As depicted in Fig. 3, a decline of 
ionic motion at even higher temperatures up to 900 K is 
indicative of Li sublattice configuration freezing, which is 
confirmed by topological structural analysis. The Arrhe¬ 
nius plot (Fig. 4) of diffusion coefficients across the tem¬ 
perature range of the CMD simulations shows a clear in¬ 
dication of a sublattice phase transition where the Li-ion 
arrangement at low temperatures freezes into its ground- 
state configuration at lower temperatures, while at high 
temperature the sublattice melts and disorder prevails. A 
relatively low activation energy Ea of 0.27 eV is obtained 
for high temperatures (1700-1100 K). But for a range of 
lower temperatures (900-500 K), the extracted Ea is as 
high as 0.98 eV. We expect the actual experimental acti¬ 
vation energy to be within these limits depending on the 
sample preparation procedure of the x = 6 crystal, that 
may result in a partially locked Li arrangement. For con¬ 
centrations just above and below x = 6, mobile defects 
move freely in the background of the stable configura¬ 
tion, leading to sharp increases of conductivity. In this 
picture the concentration of free carriers remain smaller 
that the nominal Li concentration x in this and other 
composition regions close to an ordered ground state. 
This picture explains the recent finding that only about 
10% of Li in garnets contribute to conductivity [22], and 
also suggests that half-integer concentration are optimal 
for conductivity. For x = 3 and x — 7 (t-LLZ) the CMD 
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simulations do not show any significant ionic motion even 
at high temperature. Indeed, for x = 3, an ordered dis¬ 
connected network of only tet Li cannot move because 
migration to an oct site is very unfavorable. As more 
Li ions are introduced, they act as free carriers moving 
in the background of a static high-symmetry x = 3 con¬ 
figuration, leading to a rapid increase in conductivity. 
A similar situation prevails for the fully filled set of oct 
sites in the ordered tetragonal ground state of x = 7 (t- 
LLZ), where the cost of creating oct vacancies is high. 
Reducing Li content by adjusting the composition intro¬ 
duces mobile vacancies that increase conductivity. At 
the same time, a higher degree of motion can also be ob¬ 
served in the “excited” configuration of LLZ with cubic 
structure c-LLZ and a tetioct occupancy ratio of 1.5:5.5. 
In this case oct vacancies are present and act as carriers, 
and the availability of multiple degenerate configurations 
leads to disorder and significant freedom of motion. The 
simulations predict a distinct difference of 1-2 orders of 
magnitude in the conductivity between the ground state 
(t-Li7) and the disordered state (c-Li7), consistent with 
experimental observations [23]. CMD simulations for the 
case of x = 5 exhibit better transport, as compared to 
the other integer compositions, due to a greater intrinsic 
propensity towards disorder. This can be attributed to 
the presence of multiple configurations that are energeti¬ 
cally similar to the ground-state arrangement and tend to 
yield an average low-symmetry configuration. The slight 
dip of conductivity at x = 4 is again consistent with the 
presence of an ordered ground state; however, in this case 
disorder is seen to set in at lower temperatures than for 
x = 6. 



1000/T (1/K) 

FIG. 4. Diffusion coefficient versus inverse temperature for 
the x = 6 garnet structure, illustrating the order-disorder 
phase transition at ~900K, as computed using CMD. Black 
circles represent high-temperature data and green triangles 
represent low-temperature data 

In conclusion, we provide the basis for understanding 
the non-trivial relationship between composition, struc¬ 


ture and transport that governs performance of garnets in 
the wide range of compositions. Through novel system¬ 
atic structural symmetry analysis based on the universal 
occupancy trend, combined with ab-initio energy com¬ 
putations, we predict that ordered ground state crystal 
structures exist for integer compositions. Some of these 
structures have been confirmed experimentally, while 
others are yet to be detected. In particular, we iden¬ 
tify ordered ground states at most integer compositions 
(notably new phases for the x = 4 and x = 6), and char¬ 
acterize their configurational excitations. The lack of sta¬ 
ble ordered states at half-integer concentrations results in 
disordered frustrated configurations that maximize car¬ 
rier concentration and mobility. This is a consequence 
of the strong Li-Li interaction and the particular sym¬ 
metry structure of the host lattice, and indicates that it 
is primarily the Li concentration and arrangement that 
controls transport in garnets. Based on this new under¬ 
standing, we expect that intrinsic conductivity is highest 
at compositions around x = 5.5 and x = 6.5, and that is 
where future experimental efforts should focus. We spec¬ 
ulate that presented methods and understanding are rel¬ 
evant to other families of ionic conductors, where Li con¬ 
centration is adjustable and Li-Li interaction is strong, 
such that sublattice frustration and disorder can be used 
to optimize ionic transport. 
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SUPPORTING INFORMATION: EFFECTS OF 
SUBLATTICE SYMMETRY AND FRUSTRATION 

ON IONIC TRANSPORT IN GARNET SOLID 
ELECTROLYTES 

COMPUTATIONAL METHODOLOGY 

Ab-initio simulations were performed in the framework 
of density functional theory (DFT) as implemented in the 
Quantum Espresso software [1]. Ultrasoft pseudopoten¬ 
tials [2] were used to represent ions, and the electronic 
wave functions were expanded on a plane-wave basis set 
with cut-off energy of 30 Ry. The Brillouin zone was 
sampled using only the T-point, which was sufficiently 
accurate for unit cells of about 100 atoms. The Kohn- 
Sham equations were solved self-consistently, employing 
a generalized gradient approximation [3], with a conver¬ 
gence criterion of 10' 7 Ry for the total energy. Structural 
relaxations were performed using the BFGS algorithm 
[4, 5] with the convergence threshold for forces of 10“ 3 
Ry/bohr. In ab-initio molecular dynamics (AIMD), the 
equations of motion were integrated using the leapfrog 
Verlet algorithm [6]. A time-step of 3 fs was used, and a 
Nose-Hoover thermostat [7] was applied to regulate the 
temperature range (500 - 1300 K) and to reproduce the 
thermodynamics of a canonical (NVT) ensemble. Inter¬ 
atomic interactions were modeled with rigid-ion poten¬ 
tials which consist of the Coulomb interaction and short- 
range Buckingham potentials for Li-0 and 0-0 inter¬ 
actions [8] and La-0 and Zr-0 interactions [9]. In or¬ 
der to obtain good thermodynamics and hopping statis¬ 
tics the classical molecular dynamics (CMD) simulations 
were performed in a 2 x 2 x 2 cubic supercell for a time pe¬ 
riod of 10 ns with a time-step of 1 fs. A pre-equilibration 
simulation of 2 ns was conducted prior to the 10 ns pro¬ 
duction run. At 300 - 400 K, the diffusion coefficients 
(D) extracted from these simulations are small (D < 10- 
11 cm 2 /s) and are associated with large relative errors. 

In all CMD simulations the cell geometry was taken to 
be cubic and the cation composition was fixed; only the 
concentration of Li-ions was varied in order to study its 
influence on diffusion. Due to the approximate potentials 
and non-stoichiometric compositions used, the absolute 
values are not expected to describe the actual diffusion 
coefficients quantitatively; only qualitative comparison 
between different compositions is expected to be relevant. 
In contrast, in AIMD, the cation types were explicitly 
adjusted to charge-balance the cell. 

Equation (1) shows the Einstein relation that allows 
extraction of the diffusion coefficient ( D ) from the mean- 
squared displacement of the Li + ions. The mean square 
displacement (MSD) tracks the extent of diffusion of Li + 
within the crystal by calculating the deviation of the po¬ 
sitions of the Li + at each time-step relative to a reference 
position. 


D= lim l(|r(*)-r(0)| 2 > (1) 

x^oo or 

Based on the Einstein relation, the diffusion coefficient 
is obtained from a linear fit of the MSD. The slope- 
extracted diffusion coefficient is related to the activation 
energy based on the empirical Arrhenius equation (2). 

D = D 0 e { ^ ) (2) 

The ionic conductivity (a) was calculated from the dif¬ 
fusion coefficient (D) using the Nernst-Einstein equation 
as, 

a = nq kSf (3) 

Here ks is the Boltzmann constant, T the temperature, 
n the concentration of Li + ions whose charge is q = 1. 


SYMMETRY CLASSIFICATION OF SITE 
OCCUPANCIES 

Fig. 1 depicts the generalized symmetry classification 
of the garnet crystal. In the ideal cubic structure (group 
230) Li + can potentially occupy three types of sites: oc¬ 
tahedral (6-fold coordinated with O), tetrahedral (4-fold 
coordinated), and prismatic (6-fold coordinated). Ac¬ 
cording to our DFT calculations, the structures where 
prismatic sites are occupied by Li are too high in en¬ 
ergy to be relevant, so we do not consider them in our 
analysis. For the subsequent crystallographic subgroup 
projections, a few imperative symmetry considerations 
were made: 

1. Starting with the highest symmetry, a high- 
symmetry Li + ordering is selected corresponding to 
the desired Li + content. For instance, in the case 
of Li 7 La 3 Zr 2 0 12 , the 244-atom unit cell has 56 of 
the possible 72 octahedral and tetrahedral sites oc¬ 
cupied. We then look for maximal subgroups of 
index 2 that would split the 48g and 24 d Wyckoff 
positions into inequivalent groups of smaller multi¬ 
plicities in such a way as to satisfy the occupancy 
ratio constraint (6 oct : 1 tet in the case of x=7). 
Symmetry sub-group analysis for x = 7 reveals that 
it is impossible to find a cubic subgroup for the de¬ 
sired Li + occupancy but rather the highest sym¬ 
metry attainable corresponded to the tetragonal 
Li 7 La 3 Zr 2 0 12 (t-LLZ) structure (group No. 142) 
and several smaller tetragonal subgroups. This is 
also the case for x = 5 and x = 4. 
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2. For each Li + configuration corresponding to a par¬ 
ticular subgroup and conjugacy class, a full atomic 
and crystal cell relaxation was performed using 
DFT. An energy comparison of the configurations 
allowed detection of the ground-state configuration 
(lowest-energy configuration) that is likely to be ob¬ 
served in experiment. We find that the lowest en¬ 
ergy in all the ordered configurations identified in 
this way appear in the highest subgroup of parent 
prototype group No. 230 that satisfies the occu¬ 
pancy ration constraint. The excitation energies 
are estimated as the energy difference between the 
ground state and the next lowest-energy ordered 
configuration. 



FIG. 1. Group-subgroup relationship tree that exemplifies a 
typical set of paths from the parent group 230 to lower symme¬ 
try subgroups. T and O represent tetrahedral and octahedral 
groups of sites, respectively, together with the multiplicities of 
each group. This set of relationships is not complete and only 
serves as an illustration of the symmetry projection procedure 
used to enumerate all possible ordered states. 


GROUND STATE STRUCTURES 


In Tables I II and III we list the crystallographic rep¬ 
resentations of the ground state structures for x = 4 and 
x = 6 identified through the above procedure, combining 
symmetry projections and DFT optimization. 


TABLE I. x = 9 fully filled cubic garnet prototype structure 
(stoichiometry not balanced). Group No. 230. Unit Cell: 
12.8 12.8 12.8 90.0 90.0 90.0. 

Li 24 d 0.250000 0.875000 0.000000 

Li 48 g 0.125000 0.682600 0.567400 

A (La) 24c 0.125000 0.000000 0.250000 

B (Zr) 16a 0.000000 0.000000 0.000000 

O 96/i 0.279650 0.105640 0.198940 


TABLE II. x m 6 cubic ordered state (Li 6 La 3 ZrTa0 12 ) 
Group No. 198. Unit Cell: 12.96 12.96 12.96 90.0 90.0 90.0 
Li(l)t 12 b 0.62809565 0.51777692 0.24341914 
Li(2)o 12 b 0.42134654 0.3895379 0.6889208 
Li(3)o 12 b 0.86249175 0.33228876 0.43854475 
Li(4)o 12 b 0.55829781 0.69435057 0.34323682 
La(l) 12 b 0.12753211 0.0059228 0.25388513 
La(2) 12 b 0.88209927 0.99745617 0.74645235 
Zr(l) 4a 0.0072975 0.0072975 0.0072975 
Zr(2) 4a 0.7484067 0.2484067 0.2515933 
Ta(l) 4a 0.24121324 0.74121324 0.75878676 
Ta(2) 4a 0.50141341 0.50141341 0.50141341 
0(1) 12 b 0.30114506 0.12687324 0.19004257 

0(2) 12 b 0.85095734 0.52151886 0.04205643 

0(3) 12 b 0.72512746 0.91076796 0.80135235 

0(4) 12 b 0.13773853 0.46199925 0.94492937 

0(5) 12 b 0.78156338 0.59826182 0.71760309 

0(6) 12 b 0.34127734 0.02200982 0.57727936 

0(7) 12 b 0.22591701 0.39943203 0.31618047 

0(8) 12 b 0.64875228 0.97644017 0.48229974 


TABLE III. x = 4 tetragonal ground state. Group No. 
Unit Cell: 12.77 12.77 12.64 90.0 90.0 90.0 


95. 


Li(l) 

Li(2) 

Li(3) 

Li(4) 

Li(5) 

Li(6) 

La(l) 

La(2) 

La(3) 

La(4) 

La(5) 

La(6) 

Ta 

W 

0 ( 1 ) 

0 ( 2 ) 

0(3) 

0(4) 

0(5) 

0 ( 6 ) 

0(7) 

0 ( 8 ) 

0(9) 


4a 0.000000 0.37862317 0.5000000 
4b 0.500000 0.87135575 0.0000000 
8 d 0.59084348 0.67171652 0.44378762 
4a 0.87671156 0.000000 0.7500000 

8 d 0.75902292 0.24510096 0.12652444 
4c 0.42807484 0.42807484 
4b 0.88433456 0.500000 

4a 0.62821799 0.000000 


4c 0.75830768 0.75830768 
4b 0.62810149 0.500000 

4a 0.87522166 0.000000 

4c 0.74780442 0.25219558 
8 d 0.750000 0.500000 

8 d 0.750000 


0.1250000 

0.7500000 

0.7500000 

0.6250000 

0.2500000 

0.2500000 

0.375000 

0.500000 

0.000000 


0.000000 

8 d 0.03199317 0.60993705 0.70125429 
8 d 0.47879431 0.10228318 0.8045714 
8 d 0.95223418 0.78006655 0.60957887 
8 d 0.44813813 0.71893058 0.40179069 
8 d 0.86011885 0.69897432 0.77979717 
8 d 0.64868815 0.20280054 0.72515602 
8 d 0.47251329 0.3968485 0.29770294 
8 d 0.02826262 0.89200192 0.19808783 
8 d 0.55585467 0.22400235 0.39348561 
0(10) 8 d 0.04964568 0.27821205 0.60520813 
0(11) 8 d 0.6421075 0.30037387 0.21039476 
0(12) 8 d 0.85755588 0.80240413 0.27346827 
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